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Abstract: We show that the existing CDF L = 630 pb _1 Tevatron data on pp — > ^i + ^X 
places substantial limits on a light CP-odd Higgs boson a with m a < 2m b produced via 
gg — > a, even for m a > 2m T for which BR(a — > is relatively small. Extrapolation 

of this existing CDF analysis to L = 10 fb _1 suggests that Tevatron limits on the abb 
coupling strength in the region m a > 8 GeV could be comparable to or better than limits 
from Upsilon decays in the m a < 7 GeV region. We also give rough estimates of future 
prospects at the LHC, demonstrating that early running will substantially improve limits 
on a light a (or perhaps discover a signal). In particular, outside the Upsilon peak region, 
integrated luminosity of only 5 fb _1 — 20 fb _1 (depending on m a and -y/i) could reveal a 
peak in M^+^- and will certainly place important new limits on a light a. The importance 
of such limits in the context of NMSSM Higgs discovery and {g — 2) M are outlined. 
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1. Introduction 

Many motivations for the existence of a light CP-odd Higgs boson, a, have emerged in 
a variety of contexts in recent years. Of particular interest is the m a < 2m,B region, for 
which a light Higgs, h, with SM-like WW, ZZ and fermionic couplings can have mass 
below the nominal LEP limit of nih > 114 GeV by virtue of h — )• aa —> 4r decays being 
dominant |l|, |3L [| (see also |5|, ||). For < 105 GeV, the Higgs provides perfect 
agreement with the rather compelling precision electroweak constraints, and for BR{h — > 
aa) > 0.75 also provides an explanation for the ~ 2.3a excess observed at LEP in e + e~ — > 
Zbb in the region ~ 100 GeV if rrih ~ 100 GeV. This is sometimes referred to as 
the "ideal" Higgs scenario. More generally, superstring modeling suggests the possibility 
of many light a's, at least some of which couple to r + r _ and bb. Further, it 

is not excluded that a light a with m a > 8 GeV and enhanced abb coupling could be 
responsible for the deviation of the measured muon anomalous magnetic moment a^ from 
the SM prediction |7|. Below, we will show that a light a with the required abb and a/x + ^ _ 
couplings would have been seen in existing Tevatron data for the /U + /i~ final state at low 
M^+^j-. More generally, current muon pair Tevatron data places significant limits on a 
light a. These will be further strengthened with increased Tevatron integrated luminosity 
and by data obtained at the LHC. 

The possibilities for discovery of an a and limits on the a are phrased in terms of the 
afi~ n + , ar~r + , abb and ait couplings defined via 

^/7"^//S^ 775/a - (L1) 

In this paper, we assume a Higgs model in which C aA1 - M + = C aT - T + = C^, as typified 
by a two-Higgs-doublet model (2HDM) of either type-I or type-II, or more generally if 
the lepton and down-type quark masses are generated by the same combination of Higgs 
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fields. However, one should keep in mind that there are models in which r = (C aiM -^+ = 
C ar - T +)/C abb S> 1 — such models include those in which the muon and tau masses are 
generated by different Higgs fields than the b mass. In a 2HDM of type-II and in the 
MSSM, C afl -n+ = C aT - T + = C abb = tan/3 (where tan/3 = h u /hd is the ratio of the vacuum 
expectation values for the doublets giving mass to up-type quarks vs. down-type quarks) 
and C at j = cot/3. These results are modified in the NMSSM (see, e.g. |J and ||). 1 In the 
NMSSM, both C at j and C abb = CL,-„+ = C ar - T + are multiplied by a factor cos6>a, where 
cos 6 a is defined by 

a = cos A a M ssM + sin A as, (1.2) 

where a is the lightest of the 2 CP-odd scalars in the model (sometimes labeled as ai). 
Above, cimssm is the CP-odd (doublet) scalar in the MSSM sector of the NMSSM and 
as is the additional CP-odd singlet scalar of the NMSSM. In terms of cos 6 A , C a ^—^+ — 
C aT - T + = C abb = cos 0,4 tan /3 and C at j = cos cot /3. Quite small values of cos 9 a are 
natural when m a is small as a result of being close to the U(1)r limit of the model. In 
the most general Higgs model, C a „-„+, C aT - T +, C abb and C at j will be more complicated 
functions of the vevs of the Higgs fields and the structure of the Yukawa couplings. In this 
paper, we assume C aM - M + = C ar - T + = C abb and C abb /C atJ = tan 2 /3. 

One should keep in mind, however, the fact that the above are tree-level couplings 
and that the bba^sSM coupling is especially sensitive to radiative corrections from SUSY 
particle loops that can be large when tan/3 is large 12, 13, 14]. These are typically 



characterized by the quantity which is crudely of order ■ The correction 

to the coupling then takes the form of 1/(1 + A^). Since fi can have either sign, C abb can be 
either enhanced or suppressed relative to equality with C ar - T + (the corrections to which 
are much smaller) and C a(i -u,+ (the corrections to which are negligible). 

In the past, probes of a light a have mainly relied on production of a primary particle 
(e.g. an Upsilon) which then decays to a lighter a with the emission of a known SM particle 
(e.g. a photon). Such probes are strictly limited to a maximum accessible m a by simple 
kinematics. The only exceptions to this statement have been probes based on e + e _ — > bba 
followed by a — >• t + t~ or a — > bb, with LEP providing the strongest (but still rather weak) 
limits on the a based on this type of radiative production process. 

In contrast, hadron colliders potentially have a large reach in m a as a result of the 
fact that an a can be produced via gg — > a. The gga coupling derives from quark triangle 
loops. This process, plus higher order corrections thereto, leads to a large cross section for 
the a due to the large gg "luminosity" at small gluon momentum fractions, provided the 
aqq coupling deriving from the doublet component of the a is significant. This large cross 
section will typically lead to a significant number of gg — > a — > [i^ [i~ events even though 
BR(a —7- /i + yU _ ) is not very large (and in fact is quite small for m a > 2m T ). Further, since 
the a is a very narrow resonance, all a events will typically fall into a single bin of size given 
by the M^+^- mass resolution of the experiment, typically below 100 MeV. This implies 
very controllable background levels, mostly deriving from heavy flavor production (e.g. bb 
and cc), once isolation and promptness cuts have been imposed on the muons. 



1 A convenient program for exploring the NMSSM Higgs sector is NMHDECAY Jl(], 
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The Higgs doublet component of the a can be suppressed when the a mixes with one or 
more SM singlet CP-odd field, e.g. the as of the NMSSM as made explicit above. However, 
what is important for limits is C abb . In the NMSSM context, the scenarios allowing a 
light scalar Higgs to escape LEP limits by virtue of h aa decays with m a < 2mg are 
such that cos^ is only small when tan/3 is large; for example, preferred scenarios for 
tan/3 = 10 are such that cos 9a ~ 0.1, implying C abb ~ 1. As a result, the Tevatron and 
LHC provide very significant probes of such a light a despite its being only 1% doublet (at 
the probability level). In addition, there are models beyond the MSSM [jl5| , [Tq l, including 
scenarios within the NMSSM jL7]], in which the doublet component of the a can be quite 
substantial. (For a review, see also [jf8[ .) In such models there are typically several Higgs 
scalars (including charged Higgs bosons) with masses near 100 GeV that have escaped 
discovery because of decays involving the light a. This type of scenario requires that 
tan/3 be small (1 < tan/3 < 3). When the doublet component of the a is substantial 
C abb = cos 9a tan /3 will have magnitude > 1 and hadron colliders will almost certainly 
discover or exclude the a. In the case of the NMSSM, there is a portion of the preferred 
parameter region for low tan /3 in which precisely this kind of scenario arises. However, in 
the NMSSM at low tan /3 there is a second part of the preferred parameter region in which 
there are many light Higgs bosons but the a is mainly singlet. In this latter case, \C abb \ will 
be relatively small and direct searches for the light a will be more difficult than in NMSSM 
models with larger tan/3. 

The organization of the paper is as follows. In Sec. ^, we review some basic facts about 
a light a and limits on C abb = C ail -^+ = C ar - T + coming from non-hadron-collider data. 
In Sec. H|, we discuss the additional limits that can be placed on the couplings of the a 
implied by existing Tevatron analyses and data and extrapolate these existing results to 
L = 10 fb _1 data sets. In Sec. ||[ we analyze prospects for discovering, or at least further 
improving limits on the couplings of, a light a using early LHC data. Sec. |5] summarizes 
our conclusions and provides a few additional comments. 



2. Phenomenology and limits for a light CP-odd a 

One key ingredient in understanding current limits and future prospects is the branching 
ratio for a —> [i + decays. This branching ratio (which is independent of cos 9a at tree- 
level due to the absence of tree-level a — > VV couplings and similar) is plotted in Fig. |l|. 
Note that BR{a —> /i + /x~) changes very little with increasing tan/3 at any given m a once 
tan/3 > 2. 

Limits on \C abb \ = | cos 9a | tan /3 were analyzed in (see also ^3|), based on data 
available at the time. The analysis of employed limits from T — > 7a decays, the 



importance of which was emphasized in [19| (especially within the NMSSM context), as 



well as from e + e~ — > bba production at LEP. The analysis of j7j] was done prior to the very 
recently released BaBar T(nS) — > 7a results |2^, Without including the T35 BaBar 
data, limits in the 8 GeV < m a < 2rriB range (especially, My is < m a < 2ms) are quite 
weak and suffer from uncertainty regarding % — a mixing. An update employing the ¥3$ 
data will be performed in a separate paper. In the present paper, the limits implicit in 
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Tevatron data are compared to the limits obtained in |7||. We will also briefly summarize 
how this comparison will change after inclusion of the T35 BaBar results. 

Focusing on the NMSSM, we note that it is always possible to choose cos 6 a so that 
the limits on C abb as a function of tan /3 are satisfied. The maximum allowed value of 
I cos6>a|, cos#™ ax , as a function of m a for various tan/3 values as obtained in 0] is plotted 
in Fig. ||. Constraints are strongest for m a <> 7 GeV for which Upsilon limits are strong, 
and deteriorate rapidly above that. 

Turning to the 2HDM(II), where C abb = tan/3, we note that any point for which 
cos#™ ax is smaller than 1 corresponds to an m a and tan/3 choice that is not consistent 
with the experimental limits. Disallowed regions emerge for m a <^ 10 GeV at higher tan/3. 
A disallowed region also arises over a limited m a range starting from m a > 12 GeV when 
tan (3 ~> 18, the larger the value of tan (3 the larger the interval. For example, for tan /3 = 50 
the 2HDM(II) is not consistent for m a < 10 GeV nor for 12 < m a < 37 GeV. In contrast, 
for tan /3 = 10 the 2HDM(II) model is only inconsistent for m a < 9 GeV. 

Before proceeding, we note that constraints from precision electroweak data are easily 
satisfied for a light a in both the 2HDM(II) and NMSSM cases (see @ for more discussion). 
We also wish to make note of the regions of interest for obtaining a new physics contribution, 
Aa^, of order Aa^ ~ 27.5 x 10~ 10 (the current discrepancy between observation and the 
SM prediction). These can be roughly described as follows. In the 2HDM(II) context, such 
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Figure 1: BR(a — > ) is plotted as a function of m a for a variety of tan/3 values. BR(a — > 

[i + [J.~) is independent of cos Qa at tree-level. 
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Figure 2: We plot results from for cos6^ lax in the NMSSM (where C ahl = cos 9 a tan 0) as a 
function of m a for m a > 2m T . The different curves correspond to tan ft = 1 (upper curve), 3, 10, 32 
and 50 (lowest curve). The region between ~ 2m T and ~ 8 GeV is strongly constrained by CLEO 
III data [^2| on Tig — > jt + t~ decays. The plotted limits do not include the BaBar T35 — > jt + t~ 
and T35 — » 7^i + /i _ limits that became available after the analysis of 0. 

Ja^ requires a rather precisely fixed value of tan /3 ~ 30 — 32 and m a ~ 9.9 — 12 GeV. In 
the NMSSM context, the strong constraints from Upsilon physics imply that significant 
contributions to a M are not possible until m a exceeds roughly 9.2 GeV. The maximal 5a ^ 
can exceed Aa M = 27.5 x 1(T 10 for 9.9 GeV < m a < 12 GeV if tan/3 > 32, with an almost 
precise match to this value for tan/3 = 32. For tan/3 = 50, one can match Aa^ by using a 
value of cos 9a below cos #™ ax . (The fact that matching is possible for 9.9 GeV <> m a <^ 2ms 
is particularly interesting in the context of the ideal Higgs scenario.) Further, the maximal 
<5a M is in the 7 - 20 x lO" 10 range for 12 GeV < m a < 48 GeV for tan /3 = 32 and for 
12 GeV < m a < 70 GeV for tan/3 = 50. 

At this point, it is worth discussing in more depth the "ideal" ~ 100 GeV, m a < 
2m_e, BR(h — > aa) > 0.75 Higgs scenario as discussed in [jl], ^, ||, Q. These references 
examined the degree to which obtaining the observed value of mz requires very precisely 
tuned values of the GUT scale parameters of the MSSM and NMSSM. One finds that in 
any super symmetric model this finetuning is always minimized for GUT scale parameters 
that yield a SM-like h with m/j < 100 — 105 GeV, something that is only consistent with 
LEP data if the h has unexpected decays that reduce the h — >• bb branching ratio while 
not contributing to h — > bbbb (also strongly constrained by LEP data). A Higgs sector 
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Table 1: Values of cos9a required for m a < 2ms and sufficiently large BR(h — > aa) to escape 
LEP limits on the Zbb final state. Results both without and with G < 20 required are presented 
for a selection of tan (3 values. 



tan/3 


cos 6^ ranges 


cos 9a ranges, G < 20 required 


1.7 


< -0.3 or > 0.1 


[-0.6,-0.5] or ~ +0.1 


2 


< -0.3 or > 0.1 


[-0.7,-0.5] 


3 


< -0.06 


[-0.35,-0.08] 


10 


< -0.06 or > 0.035 


[-0.12,-0.08] or [0.06,0.08] 


50 


< -0.04 or > 0.04 


[-0.06, -0.04] or ~ +0.04 



with a light a for which BR(h — > aa) > 0.75 and with m a small enough that a decays to 
BB final states are disallowed (i.e. m a < 10.56 GeV) provides a very natural possibility 
for allowing minimal finetuning. The NMSSM provides one possible example. As a useful 
benchmark, in the context of the NMSSM the tan /3 = 10 scenarios that yield the required 
m a < 2rriB and BR(h — > aa) > 0.75 are ones with 0.35 > \C abb \ (| cos 9a\ ^ 0.035). The 
lower limit arises from the fact that BR(h -4 aa) falls below the 0.75 level needed for the 
ideal Higgs scenario if | cos 6a\ is too small. Prom Fig. ^ we see that such | cos 6a\ values are 
not yet excluded for any m a > 2m T . This range becomes more restricted if, in addition, 
one requires small finetuning of the A\ and A K soft-SUSY-breaking NMSSM parameters 
that determine the properties of the a — such finetuning is characterized by a parameter 
we call G, defined in §. At tan/3 = 10, 0.6 <, \C abl \ <, 1.2 (0.06 <> \ cos6 A \ £ 0.12) is 
required if G < 20 is imposed as well as requiring m a < 2ms and BR(h — > aa) > 0.75. 
For tan /3 < 2, the means for escaping the LEP constraints on the light scalar Higgses 
are a bit more complex since two ~ 100 GeV Higgses can share the ZZ-Higgs couplings 
squared, but there is always a lower limit on | cos^l for which such escape is possible. In 
Table p], we tabulate more precisely the values of cos 9a for various tan/3 values that: a) 
have m a < 2rriB and large enough BR{h — > aa) to escape LEP limits on the h, with no 
constraint on G; and b) have small A K , A\ finetuning measure G < 20 as well as m a < 2mg 
and large enough BR(h — > aa). 

We can summarize the implications of this table as follows. First, comparing to the 
existing limits on | cos 9a \ as plotted in Fig. we see that only ideal Higgs scenarios 
(i.e. ones with rrih < 105 GeV and BR(h — > aa) large enough to escape LEP limits) 
with tan /3 > 30 and m a < 8 GeV are excluded. Ideal Higgs scenarios with tan /3 < 10 
are fairly far from being excluded. If we wish to eliminate ideal Higgs scenarios then: 
for 1.7 < tan/3 < 2, 2 we must exclude |C afe ^| > 0.17; for tan/3 = 3, we must exclude 
\C abb \ j> 0.18; for tan/3 = 10, we must exclude |C afe ^| J> 0.35 and for tan/3 = 50, we must 
exclude |C a ^| J> 2. If we only wish to exclude such scenarios that also have G < 20, then the 
required |C a6 ^| levels for tan/3 = 1.7,2,3,10,50 are 0.17,1,0.24,0.6,2, respectively. As we 

2 Scenarios with tan/3 much below 1.7 are problematical since it is difficult to retain perturbativity for 
Yukawa couplings all the way up to the unification scale. 
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shall see, completely probing even the latter levels for all m a < 2m b will be challenging, but 
hadron colliders may ultimately play a leading role. Indeed, those scenarios with G < 20 
typically have m a values above 7.5 GeV and most often above Mx 3S . Of course, the many 
scenarios with larger | cos 9a | than the values listed above will be correspondingly easier to 
exclude or verify. 

Finally, we comment on the implications of the recent ALEPH results [^4|] which place 



as a function of and m ai . 3 (Above, we use the more precise notation h\ and a\ 
for the lightest CP-even and CP-odd Higgs bosons of the NMSSM.) According to the 
ALEPH analysis, to have m hl < 100 GeV, £f < 0.52 (0.42) is required if m ai ~ 10 GeV 
(4 GeV). These limits rise rapidly with increasing m/ ll — for m/ ll = 105 GeV (the rough 
upper limit on m/ ll such that electroweak finetuning remains quite small and precision 
electroweak constraints are fully satisfied) the ALEPH analysis requires £ 2 ^ 0.85 (<^ 0.7) 
at m ai ~ 10 GeV (4 GeV). These limits are such that the easily viable NMSSM scenarios 
are ones: i) with m ai below but fairly close to 2ms, which is, in any case, strongly preferred 
by minimizing the light-ai finetuning measure G; and/or ii) with tan/3 relatively small 
(< 2). These are also the scenarios for which Upsilon constraints are either weak or absent. 
Details will be provided in a forthcoming paper [^5|]. Here, we present a simple summary. 
In particular, we note the following: a) all tan/3 < 2 cases provide m/ ll < 100 GeV 
scenarios that escape the ALEPH limits; b) there are a few G < 20, tan /3 = 3 scenarios 
with m/j 1 as large as 98 GeV and 99 GeV and with £ 2 essentially equal to the ALEPH 
limits of £ 2 < 0.42 and £ 2 < 0.45 applicable at these respective m^ values; c) tan/3 = 10 
ideal scenarios easily allow for ~ 100 — 105 GeV (because the tree-level Higgs mass 
is larger at tan/3 = 10 than at tan/3 = 3) and at m ai < 2ms many m/ ll > 100 GeV 
points have £ 2 < 0.5 in the fixed-/U scan and a few of the full-scan points have £ 2 < 0.6 for 
m/ ll ~ 105 GeV, both of which are below the m ai = 10 GeV ALEPH upper limits on £ 2 
of 0.52 at m/jj ~ 100 GeV and 0.85 at = 105 GeV; d) at tan (3 = 50 there are some 
G < 20 points with m hl ~ 100 GeV and m 0l < 2m B having £ 2 below the 0.52 ALEPH 
limit. Finally, we note that for the entire range of Higgs masses studied the ALEPH limits 
were actually ~ 2a stronger than expected. Thus, it is not completely unreasonable to 
consider the possibility that the weaker expected limits should be employed. These weaker 
limits for example allow £ 2 as large as 0.52 at ~ 95 GeV and 0.9 for ~ 100 GeV. 
These weaker limits allow ample room for the majority of the m ai < 2ms ideal Higgs 
scenarios. 

3. The role of the Tevatron 

Potentially, the Tevatron can probe precisely the m a range close to and above the T(nS) 
masses which cannot be probed in T(nS) decays. Some relevant analyzes have been per- 

3 In this arXiv version of this paper, we have access to final results and so we have updated this paragraph 
relative to the published version of this paper. 
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Figure 3: The total cross section for a production at the Tevatron is plotted vs m a for tan/3 = 
1,2,3, 10 (lowest to highest point sets). For each m a and tan/3 value, the lower (higher) point is 
the cross section without (with) resolvable parton final state contributions. 



formed looking for a very narrow resonance, denoted e, that is produced in the same way 
as the Tig. The published/preprinted results are those of [^] and |^7j from the CDF 
experiment. The latter results employ data corresponding to L = 630 pb _1 and exclude 
the potential e peak at M^+^- ~ 7.2 GeV present in the L = 110 pb _1 data of the first 
paper. The analysis was only performed for the region from 6.3 GeV < M„+„- < 9 GeV. 
The reason for not performing the analysis at lower M„+„- is that the acceptance of the 
pair relative to that of the Tig (used as a normalizing cross section) would be highly 
mass dependent. This is due to the fact that CDF is only able to see pairs with 

Pt > 5 GeV and for M^+^- < 6.3 GeV the fraction of pairs that fail this cut becomes 
highly mass dependent. No reason for not analyzing the region of M^+^- > 9 GeV is given, 
although it is in this region that the T(1S,2S,3S) peaks are present. 

Our goal here is to use the above e analyses to place limits on a CP-odd a. This is 
possible under certain assumptions detailed below. In particular, we will place limits on 

^abb- 

The dominant production mechanism for a light a at a hadron collider is different 
than that for the e (assumed to be the same or very similar in kinematic shape etc. for 
the Tig). The production mechanisms for the Tig remain uncertain. It has recently been 
claimed p8| that an NNLO version of the leading order (LO) calculation can reproduce 
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the Tevatron results for direct production of the T15 at larger px (roughly p? > 5 GeV). 
The diagrams employed begin with the LO O(a^) process gg — > bbg with the bb pair 
turning into the T15 with probability determined by | i?x is (0) | 2 , leading to an T\s + g final 
state. At NLO, the a% diagrams include virtual correction diagrams that also lead to the 
YiS + 9 final state and several diagrams containing an extra quark or gluon in the final 
state (T15 + 2g and T15 + bb). Several C(a|) diagrams leading to T15 + 3j final states 
(especially T15 + 3g) are argued to be of importance at larger pt and are also included. 
Resummation [^9| is necessary to get the low px portion of the cross section. From CDF 
and DO data, the direct Tig production cross section is measured to be about 50% of the 
total. Indirect contributions coming from, for example, gg — > \b followed by \b TiSl 
make up the remaining 50% of the total T15 production rate. In contrast, a, being a 
spin-0 resonance, will be dominantly produced via gg — )• a through the quarkdoop induced 
gga coupling. In addition, there are large QCD corrections to the onedoop-induced cross 
section. These are of two basic types: a) virtual corrections and soft gluon corrections; b) 
corrections containing an extra resolvable gluon or quark in the final state (the dominant 
diagram is gg — > ag) in close proximity to the a. The total cross sections predicted by 
HIGLU [3C] are plotted as a function of m a for C . j = 1/C at j = tan/3 = 1, 2, 3, 10 in Fig. || 
with and without the resolvable parton final state QCD corrections. The HIGLU results 
agree well with a private program for this process. We note that the cross sections do not 
scale precisely as tan 2 /3 at large tan j3 (as naively predicted by dominance of the 6-quark 
loop diagram for the gg — >■ a coupling at high tan j3) due to the virtual corrections. In any 
case, very substantial cross sections are predicted. In the NMSSM context, at any given 
tan/3 value one should reduce the plotted result for that tan/3 by a factor of {cos 9 a) 2 . 

Among the cuts employed in the CDF analysis there is an isolation requirement 
whereby events are only included if both muons have less than 4 GeV scalar summed 
px in a cone of size AR = 0.4 about the muon. The impact of the isolation requirement 
was studied for the T15 and it was found that this isolation requirement was 99.8% effi- 
cient for the Tis despite the fact that Tig's are produced along with one or more extra 
particles in the final state. Thus, in our analysis for the a we will make the assumption 
that the components of the a cross section coming from final states containing an extra 
q or g are not significantly affected by the isolation cut and thus we will employ the full 
QCD-corrected a cross section. In addition, in the analysis of [27j only events for which 
the /i 4 ' pT pair resides in the \y\ < 1 region are retained. Thus, what we actually employ 
are the cross sections 

da(a) 

a { a )\y\<l ~ 



x2, (3.1) 

y=o 



dy 

which is an excellent approximation given that the cross section is essentially flat in y over 
this region. At the Tevatron, the ratio 



dcr(a) 
dy 



y=0 (3.2) 



o-(a)tot 

varies from roughly 0.12 at m a ~ 2 GeV to 0.19 at m a ~ 12 GeV with very weak dependence 
on tan/3. At m a = Mx is the ratio is ~ 0.15. 
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In [27 1, what is given are limits on the ratio for production of a very narrow resonance, 
the e, relative to that for the Tig 



R 



a(e)BR(e -> 
a(T ls )BR(T ls => n+n- 



(3.3) 



under the assumption that the same mechanism is responsible for e production as is respon- 
sible for Yis production. As stated earlier, since the a can be produced directly via gg — ^ a 
whereas the Tig cannot, an interpretation for the a of the limits given for the generic e 
requires actually knowing what the Tig cross section is. It also requires an assumption 
regarding the efficiency for the a of acceptance and isolation requirements relative to those 
employed for the e. 

Our analysis is the following. In Ref. [26], it is stated that the cross section for Tig 
production in the \y\ < 0.6 region at y/s = 1.8 TeV was measured to be 34,600 pb. In 
contrast, the cuts of Ref. [26 1 accept Tig events with \y\ < 1. In Ref. 26], it is stated that 
the efficiency for Tig detection (due to geometric and kinematic acceptance cuts as well 
as trigger and reconstruction efficiencies, but before imposing the isolation requirement 
noted above) is 0.066. From [31], we infer that this acceptance times efficiency factor is 
one that applies at any fixed value of y that is relatively central and after integrating over 
accepted pr's. Then, using BR(Tis — > /i + /x~) = 0.0248 and the integrated luminosity of 
L = 110 pb -1 one then predicts 

34600 x f-^j x 110 pb" 1 x 0.0248 x 0.066 = 10383 (3.4) 

events where the parenthetical fraction corrects for the increased Ay acceptance compared 
to that used in measuring the Tig cross section. This compares favorably to the 9838 
number of events that were observed before including the isolation cuts and promptness cuts 



of Table 1 in Ref. 



A cross check on the cross section is to note that the da( -J ls ^ 

ay 



BR(T 1S 
on Eq. (| 

da(Ti S ) 



fjTfJl 
of 



y=o 



753 pb value measured in [31] is comparable to the estimate based 



dy 



x BR(T 



is 



ff(Tis) 



j/|<0.6 



34600 pb 



Ay = 1.2 



x 0.0248 ~ 715 pb. (3.5) 



Because the earlier paper [31] may have employed slightly different procedures, efficiencies 
and so forth, we use the value of Eq. (3J3) at yfs = 1.8 TeV. 

Moving to the higher energy of y/s = 1.96 TeV, it is stated in [^7]] that the |y| < 0.6 
Tig cross section increases relative to yfs = 1.8 TeV by about 10%, implying 



da{T 



is) 



dy 



(1.96 TeV) x BR{T 



is 



fjt'H 



787 pb . 



(3.6) 



y=0 



This is the value we shall employ. As another cross check, we note that a 10% increase in the 
total cross section would yield about 38, 330 pb at -^/i = 1.96 TeV. At = 1.96 TeV, p7| 
states that 52,700 Tig — > events are observed using the same cuts as in Ref. |2(| (that 
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imply that only an T15 with \y\ < 1 will be accepted) and after imposing the isolation and 



promptness criteria detailed in Ref. [26|. The latter imply an additional efficiency factor of 
0.921 relative to the 0.066 efficiency referenced earlier. Multiplying these two efficiencies 
yields a net efficiency of ~ 0.061. With the 10% cross section increase and accounting for 
the increased luminosity of L = 630 pb _1 , the 0.061 net efficiency implies an expected 
event number of 60,244. Although this is not in perfect agreement with the 52,700 events 
actually observed, we will use the result of Eq. (3.6) below. 

Relative to the T15 efficiency, purely geometric effects alter the efficiency for e pro- 
duction and we assume that the same geometric changes apply to the a. The formula of 
|H is: 



efficiency(e) = ef ficiency(T ls ) 



0.655 + 



(0.974 - 0.655) 
(9.0 - 6.3) 



(m € — 6.3) 



(3.7) 



We will employ this same relative efficiency for the a as a function of m a using as well 
ef ficiency(Yis) = 0.061 as obtained above. 

The most precise limits on C abb are obtained using the ratio R defined in Eq. 
We recall from Refs. I 



27] that the limits on R are obtained by performing a smooth fit 



to the event distribution and looking for fluctuations about this smooth fit. The limits on 
the e (or the a in our case) are then obtained by placing small Gaussians at each possible 
m a value and placing limits using the observed fluctuations about the smooth fit. In the 
mass region for which CDF has performed this analysis, 6.3 GeV <> m a <> 9 GeV, it is 
very convenient to simply directly employ their results. As stated, we assume that the a 
efficiencies are the same as for the e, in which case we can compute the ratio R as 



R ~ 



dcr(a) 
dy 



y=0 



x BR(a — > n + fjL ) 



dy 



x BR(T 1S -> m+aO 



(3- 



where 

dv(T ls ) 
dy 



dcr(a) 
dy 



BR(T IS 



is computed using HIGLU, BR(a 
is as given in Eq. 



) is taken from Fig. |l] and 
Note that the exact values of the 



y=0 

efficiencies, Eq. fl3.?1), are not important using this procedure so long as the efficiency for 
the a is the same as for the e. 

With these assumptions and inputs we can then predict the ratio R for the case of 



e = a and compare to the 90% CL upper limits of [[F?]] based on L = 630 pb^ 1 of analyzed 
CDF data. This comparison appears in Fig. |] for a number of C , r = cos 9a tan (5 choices. 
We observe that the predicted R depends almost entirely on |C a6 ^|, with extremely little 
dependence on tan/3 separately for the tan/3 > 1, m a > 4 GeV parameter region on 
which we focus. The corresponding bin-by-bin limits on |C afe ^| obtained by interpolation 
appear in Fig. §. In the 2HDM(II), they are limits on tan/3 = C ahl . In the NMSSM, 
these are limits on C af r = cos 9a tan (3. In both cases, the interpolations are only accurate 
for tan /3 > 1 and m a > 4 GeV. From Fig. ||, we find that the limits based on the 
existing L = 630 pb _1 analysis roughly exclude \C b r\ > 3 for 6.8 < m a < 9 GeV and 



C abl \ > 2 for 8.2 < m a < 9 GeV, but do not exclude C, 



abb 



1 for any of the m a values 



- 11 - 



Tevatron Di— muons 




2 4 6 8 10 12 

m a (GeV) 

Figure 4: We plot the 90% CL limits on the ratio R for a production at the Tevatron as a function 
of m a compared to NMSSM predictions using HIGLU for the following cases: (tan/3 = 1,cos9a = 
1) (red +'s), (tan/3 = 2 1 cos9a = 1) (blue diamonds), (tan/3 = 3,cos#a = 1) (green x's) and 
(tan/3 = IOjCOS^ = 0.1) (yellow squares). 

in the analysis range. In the 2HDM(II) case the C abb = tan /3 limits from the Tevatron 
are stronger than those from Upsilon decays and LEP data, as summarized earlier, for 
8 GeV < m a < 9 GeV. In the NMSSM case the limits on |cos#a| from the Tevatron 
data are the stronger in much the same mass range, as we detail shortly. In Fig. |5| we 
also plot the statistically extrapolated limits that would result by increasing the data 
sample to L = 10 fb _1 . (Presumably, a real analysis of a high luminosity data set would 
do even better.) Since the a signal cross section varies roughly as (C abb ) 2 , even this large 
luminosity increase leads to limits that are improved by only a factor of a bit more than two. 
Nonetheless, one approaches the \C abb \ ~ 1 level of interest in the NMSSM at m a = 9 GeV. 

Focusing now on the NMSSM, we compute the upper limit on cos 6 a, cos #™ ax , obtained 
by appropriate interpolation of the results of Fig. |I[ We again emphasize that although 
BR(a — > /x + /i _ ) is tan/3 dependent as shown in Fig. [l], it is nonetheless the case that 
for tan (3 > 1 and m a > 4 GeV the limits on C abb = cos 6 a tan f3 are almost independent 
of tan/3 at fixed C abb , as found in Fig. || (compare the two cos 6a tan /3 = 1 cases - 
tan/3 = cos #,4 = 1 vs. tan/3 = 10, cos9a = 0.1). This can be understood as follows. At 
low tan/3, although BR(a — > fi + is suppressed, contributions to the gga coupling from 
loops involving the top quark are substantial relative to loops involving the bottom quark. 
In comparison, at large tan/3 one finds that BR{a — > is maximal but top-quark loops 

are relatively suppressed compared to bottom quark loops. These two effects very nearly 
cancel one another leaving the net a cross section unchanged at fixed C bb = cos 6a tan /3. 
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Figure 5: We plot the 90% CL upper limits on \C ab j\ obtained using the results for the ratio R 
of Fig. g. The 10 fb -1 results are obtained by statistical extrapolation of the 630 fb 1 results. In 
the context of the 2HDM(II), C abl = tan/3. In the context of the NMSSM, C abl = cos 6 A tan /3. In 
both cases, limits were derived assuming tan/? > 1. 



As a result, it is easy to extract the cos#™ ax values for different tan/3 values directly from 
the plotted limit on \C abb \ shown in Fig. [|. 

The resulting cos#™ ax limits are shown in Fig. ^ in comparison to the upper limits 
plotted earlier in Fig. ^. The figure focuses on the 6 GeV <^ m a <^ 9 GeV region for which 
we have extracted the Tevatron limits using R. What we observe is that the 630 pb -1 90% 
CL limits become the strongest for m a > 8.3 GeV. In a forthcoming paper, we will analyze 
the impact of BaBar data for T35 — > jt + t~ decays on the 6 GeV < m a <^ 10 GeV region. 
Our preliminary results suggest that the Tevatron limits plotted above and the T35 limits 
are very similar at m a ~ 9 GeV, with T35 limits being superior for lower m a . 

Given the above, the great value of extending the Tevatron analysis above M M +^- = 
9 GeV is apparent. A full analysis of existing and future data all the way out to m a ~ 
12 GeV is needed. First, it might strongly constrain the properties of any light a with 
m a < 2rriB that would allow for the ideal Higgs scenario. Second, it might completely 
eliminate the possibility that a light a could provide a major contribution to a^. At the 
moment, the Tevatron limits on C abb shown rule out a significant contribution to Aa M from 
an a with m a < 9 GeV, while in the range 9 GeV < m a < 2m b these limits leave open the 
possibility of Ao^ arising from diagrams involving a CP-odd a. Only the Tevatron and/or 
LHC can probe the region of m a above the Upsilon masses. 

Absent a full analysis by CDF of limits on R in the region M^+^~ > 9 GeV and given 
that this region is of great interest, we wish to make some estimates of limits on \C 
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Figure 6: We plot the 90% CL cos 6^ values in the NMSSM context obtained from the results 
of Fig. |i| for the 630 pb _1 CDF data set, in comparison to the cos 6*™ ax values plotted in Fig. 0. For 
clarity, the plot is limited to the m a region over which the Tevatron data are relevant. The curve 
types are as in Fig. |[ The Fig. || results for a given curve type are those for which cos 0™ ax starts 
at lower values at low m a rising to higher values at higher m a . The new CDF limits are those that 
begin near m a ~ 6.3 GeV and terminate at m a ~ 9 GeV and that fall (with fluctuations) as m a 
increases. 



based on the event number plots of Ref. [27]. We have employed the following procedure. 



First, for this analysis, we must know the efficiency for detecting the a. For our estimates 
we use ef ficiency(m a ) from Eq. Q3.7|) and ef ficiency{T \s) = 0.061 (as motivated earlier 
below Eq. ( |3.6| )) to predict the number of a events as a function of m a - Second, we wish 
to determine how many of the total number of a events fall into a 50 MeV bin centered 



on m a . To do so, we need to know the resolution as a function of m a . In [27|, it is stated 
that the resolution, a r varies from 32 MeV to 50 MeV in going from m a = 6.3 GeV to 
m a = 9 GeV, with a value of 52 MeV at Mx is . We use a simple linear interpolation for 
other values of m a , but do not allow a T to fall below 25 MeV at low m a . The fraction of a 
events distributed as a Gaussian of width a r that fall into a 50 MeV bin (which should be 
thought of as a bin of half width 25 MeV) that is centered on m a is given by 

f(m) = Erf (^MfL) (3.9, 

where a r (m) is in MeV. In Fig. [?], we plot the 1.646cr, i.e. 90% CL, fluctuation number for 
each of the CDF 50 MeV bins compared to the predicted number of a events that would 
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Figure 7: We plot the upwards fluctuation in the number of events in a given bin corresponding to 
a 1.646cr (90% CL) excess as predicted using the L — 630 pb -1 event numbers of p7| . These limits 
are compared to the predicted number of events for an a resonance centered on that bin spread out 
by the experimental resolution. The same (tan/3, cos 9 a) values as in Fig. |] are considered. 



fall into that bin. We do this for the same selection of (tan/3, cos 9a) values as employed in 
Fig. ^| One observes that for m a ~ 6 GeV (m a ~ 9 GeV) 90% CL sensitivity is anticipated 
for C abb J> 3 {C abb J> 2). This anticipates in an average sense the more precise (and more 
fluctuating) results based on the R analysis found in Fig. |5[ 

Sensitivity to the a in the S/^/B sense can actually be improved by taking a bin 
size that properly matches a r . If the background is flat then the optimal bin size is 
which retains a fraction Erf (I) = 0.843 of the total a signal and yields B = 2\^2a r d ^" B — • 
Following this procedure we can then use interpolation to extract the \ C abb \ value such that 
S/\^B is 1.646. The resulting values of \C abb \ which correspond to this 90 CL fluctuation 
in the AM At + M - = 2\f2a r acceptance window are plotted in Fig. ||. As anticipated above, 
this event counting method turns out to give a good average representation of the results 
obtained using R (which analysis was based on bin-by-bin fits of the fluctuations about 
a smooth curve) at the 90% CL. Thus, despite relatively small S/B levels (typically of 
order 0.02 in each of two neighboring bins for a 1.646cr net fluctuation), our estimates 
for expectations for m a > 9 GeV (using the approach of assuming there were no 1.646cr 
fluctuations in the absolute L = 630 pb -1 event numbers in acceptance windows of size 
AM^+^j-) should give a good idea of the limits that are implicit in current data. 
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Figure 8: We plot approximate limits on \C b r\ as a function of m a estimated assuming that 
bins centered on m a and encompassing an m a range of 2^2 x ay, where oy is the experimental 
-^i+/x- resolution at the given m a , do not have a 90% CL fluctuation relative to the number of 
events observed and plotted over this bin range in M for L = 630 pb . The L = 10 fb _i 
histogram corresponds to simply scaling the predicted a event rate and the 90% CL fluctuations to 
the increased luminosity. 



As an aside, we note that even though the shape of the Y(nS) resonances (which are 
also very narrow) will also be determined by oy, one can learn if there is an excess in 
the p + fj,~ final state by also looking at the T(nS') resonance in the e + e~ final state to 
which the a will not contribute. Assuming lepton universality, the T{nS) contribution in 
the n + pT final state can be subtracted from the p + p~ spectrum, after which any residual 
excess from the presence of an a would become apparent. Statistical errors resulting from 
this subtraction will be roughly a factor of \/2 larger than employed above. However, this 
procedure does rely on a precise understanding of efficiencies, resolutions and the like for 
electrons. Another technique that could be considered is comparing one Upsilon resonance 
to another. If the relative normalization between the two Upsilon resonances can be suffi- 
ciently precisely predicted, including both theoretical and experimental uncertainties, then 
an a —> p + p~ signal hiding under one of the Upsilons could become apparent. 

Both CDF and DO will continue to accumulate data far in excess of L = 630 pb -1 . 
Thus, it is useful to extrapolate to higher luminosity using the observed number of events 
for L = 630 pb -1 plotted in [27|. We rescale the observed number of events in each 
bin to L = 10 fb _1 and compute the 90% CL fluctuation upper limit in each bin as 
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Figure 9: We plot the upwards fluctuation in the number of events in a given bin corresponding to 
a 1.646c7 (90% CL) excess as predicted using the L — 630 pb^ 1 event numbers of |27| scaled up to 
L = 10 fb -1 . These limits are compared to the predicted number of events for L — 10 fb -1 for an a 
resonance centered on that bin spread out by the experimental resolution. The same (tan j3, cos 6 a) 
values as in Fig. || are considered. 

1.646 x 1 N ev t(bin). In Fig. ^, we plot these extrapolated 1.646cr fluctuation levels in each 
bin compared to the predicted number of a events in each bin for the same selection of 
(tan /3, cos 9 a) values as employed in Fig. The extracted limits on \C ab ^\ are plotted in 
Fig. ||. We see that the eventual Tevatron limits from just one experiment could easily 
be superior to those currently available from Upsilon decays, in particular probing the 
C a bb = cos ® A * an P ~ 1 coupling level, of particular interest for the "ideal" Higgs scenario, 
all the way up to m a = 2m#, except in the vicinity of the Upsilon peaks. 

Some further comments are the following. First, we emphasize that the 1.646<r ap- 
proximate procedure leads to the expectation that quite important limits are possible for 
L = 10 fb -1 in the 9 GeV < m a < 2mg region of great interest in the NMSSM version 
of the ideal Higgs scenario for which m a > 8 and |C afe ^| ~ 0.2 — 2 is a strongly preferred 
parameter region. If this scenario were to be nature's choice, there is a decent chance of 
observing an o using the dimuon spectrum analysis and the ultimate Tevatron data set. 

Second, we again note that even the L = 630 pb _1 estimated limits of Fig. |8| for 
9 GeV < m a < 12 GeV would rule out the enhanced \C ab ^\ values of order 30 needed for 
a-exchange graphs to explain the discrepancy for m a in this mass region, which is the 
only relatively low mass region for which other current constraints are sufficiently weak 



-17- 



that the a might provide the observed discrepancy. It is thus quite important for CDF 
(and DO) to perform the needed analysis using the R or similar technique and determine 
whether or not the rough limits we obtained above using event numbers are approximately 
correct. 



4. LHC Prospects 

The basic question is whether the LHC will be able to improve over the Tevatron L = 
10 fb _1 projected results and limits. The total cross sections at the LHC appear in Figs. 10 



and 11 for ^/s = 14 TeV and y/s = 10 TeV respectively; they are plotted analogously to 
those for the Tevatron appearing in Fig. ^. Recall that these cross sections are those 
appropriate in the 2HDM(II) context. We see that, relative to the Tevatron, the \fs = 
14 TeV cross sections are about a factor of 3 — 7 higher, the smaller (larger) ratio applying 
at small (large) m a . Relative to the -y/i = 14 TeV cross section, the -^/i = 10 TeV cross 
sections are roughly a factor of ~ 1.2 smaller at m a = 2 GeV and a factor ~ 1.34 smaller 
at m a = 10 GeV, more or less independent of the tan/3 value. It now appears that perhaps 
as much as a year will be spent running at -^/s = 7 TeV. Thus, we also plot in Fig. [l^ the 
cross sections for this latter energy. At m a = 2 GeV, the a cross section at \/s = 10 TeV is 
a factor of about 1.15 larger than the yfs = 7 TeV cross section; the factor rises to ~ 1.35 
for m a = 10 GeV. The modest decrease of the a cross section with decreasing energy is a 
result of the fact that the gg luminosity at low m a varies slowly with y/s. This is one of 
the reasons why searches for a light a are very appropriate in early LHC running. 

In going to the NMSSM, one takes these results for any given tan f3 choice and then 
multiplies by (cos 9a) 2 , the square of the overlap fraction of the a with the 2HDM compo- 
nent. 

As an example of how limits obtained at the LHC in early running will compare to 
the Tevatron limits, let us consider the case of tan/3 = 10 and cos^ = 0.1, for which 
C ab j; = 1. As shown in Fig. B even with L = 10 fb _1 the Tevatron is not fully able to 
probe at the 90% CL the predicted relatively small a event levels except at m a values close 
to 2rriB but outside the T(nS) peaks. In more detail, for the above parameter choices, the 
predicted number of a — > events for \y\ < 1 in a AM^+^~ = 2\[2o r bin centered on 

m a is 436, 615 and 475 at m a = 8 GeV, Mx is and 10.5 GeV, respectively, where the event 
numbers quoted incorporate the Erf{l) = 0.8427 reduction factor associated with keeping 
only events in an interval of size AM^ +At - = 2y/2a r . The actual AM^+^- = 2^2a r 
values are 43 MeV, 52 MeV and 57 MeV at 8 GeV, M Tls and 10.5 GeV, respectively. 
As regards the background, we take the 50 MeV bin event numbers in the CDF plot of 
the number of events in each bin and rescale to the AM^^- interval sizes at the above 
Mp+fj,- = m a choices. This gives us a background event number A^AAf fi+fi _ at each m a . 
The la fluctuations in these background event numbers, . N&m + _ , are 468, 945 and 285, 
respectively. The statistical significances of the a signals are then ~ 0.93cr, ~ 0.65<r and 
~ 1.67cr, respectively. Only the latter is (slightly) above the 1.646cr level corresponding 
to 90% CL. However, to repeat, this high m a < 2rriB region is particularly favored in the 
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Figure 10: The total cross section for a production at the LHC for y/s = 14 TeV is plotted vs m a 
for tan j3 = 1, 2, 3, 10 (lowest to highest point sets). For each m a and tan /3 value, the lower (higher) 
point is the cross section without (with) resolvable parton final state contributions. 



model context. But, to reach 5a at m a = 10.5 GeV would require about nine times as much 
integrated luminosity, i.e. L ~ 90 fb^ 1 and 5a at m a = Mx is would require L ~ 590 fb _1 . 



Projections for the LHC have been made public by ATLAS. In Fig. 1 of [32], one 
finds a plot of da/dM fl + fl - coming from bb production, Drell-Yan production and Tig 
production. The dimuon Drell-Yan contribution is negligible compared to that from bb 
production even after the latter is reduced by muon isolation requirements. We ignore the 
Drell-Yan contribution in all subsequent discussions. 

In generating the bb and Tig cross sections, only events with px cuts requiring one 
muon with px > 6 GeV and a 2nd muon with > 4 GeV, both with \r]\ < 2.4, were 
retained. A recent Monte Carlo study [33] finds that these events constitute 20% of the total 
inclusive cross section. The fraction of these events that survive after further requirements 
related to triggering, reconstruction and the final analysis selection cuts is 50%. Thus, the 
net efficiency for the T\s events plotted in Fig. 1 of (3^] is ~ 0.5 x 0.2 = 0.1. Therefore, we 
will write cat las = O.lr for the fraction of inclusive a events that will be retained, where 
r ~ 1 for the cuts and triggering strategies studied so far, but r > 1 is probably achievable 
if these are optimized for the CP-odd a. 

Returning to Fig. 1 of |32|, we observe a 66-induced dimuon cross section level for 
y/s = 14 TeV of order da/dM^^- ~ 50 - 90 pb/100 MeV in the M^+^- £ [8 GeV, 2m B ] 
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Figure 11: The total cross section for a production at the LHC for y/s = 10 TeV is plotted vs m a 
for tan j3 = 1, 2, 3, 10 (lowest to highest point sets). For each m a and tan /? value, the lower (higher) 
point is the cross section without (with) resolvable parton final state contributions. 



interval when outside the Upsilon peak region. This is the dimuon cross section from bb 
heavy flavor production only. The author of [32] estimates |34|] that one should at most 
double this cross section to account for cc production and other contributions. We will 
make estimates based on multiplying the 66-induced dimuon cross section by a factor of 
two. To this, we add the T\s cross section as plotted in Fig. 1. The net resulting spectrum 
constitutes the background to the a signal that we discuss shortly. 

As in the CDF case, we will use a bin size of Am + _ = 2\[2a r (which optimizes S/y/~B 
for a flat background) for comparing the a signal to the above stipulated background. As 
for resolutions, it is stated in |^] that the resolution at the J/ip is around 54 MeV while 
that at the T±$ is close to 170 MeV. We use a linear interpolation for other values of 
M^+^-. Assuming L = 10 pb _1 of integrated luminosity, the background event numbers 

N&m H m the intervals of size AM^^- = 2y/2a T are 4055 at m a = 8 GeV, 50968 at 

m a = My is and 9620 at m a = 10.5 GeV. We take the square root to determine the la 
fluctuation level. 

We now consider the a — > /i 4 ' \x~ signal rates. From Fig. |l(], we see that at tan (3 = 10 
the total a cross section ranges from about 4.2 x 10 5 pb(cos 6a) 2 ~ 4200 pb at m a = 8 GeV 
to ~ 8500 pb at m a < 2m b for y/s = 14 TeV. The cross section for a — > fi + assuming 
tan/3 = 10 and cos6U = 0.1 will then range from 4200 - 8500 pb x (BR(a -)■ ~ 
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Figure 12: The total cross section for a production at the LHC for y/s = 7 TcV is plotted vs m a 
for tan j3 = 1, 2, 3, 10 (lowest to highest point sets). For each m a and tan /? value, the lower (higher) 
point is the cross section without (with) resolvable parton final state contributions. 



0.003) ~ 12 — 25 pb. As discussed above, we will write the total a efficiency in the form 
zatlas = 0.1 x r. Multiplying the above cross section by €atlas and by the Erf (I) = 
0.8427 acceptance factor for the ideal interval being employed and using L = 10 pb -1 
(as employed above in computing the number of background events), we obtain o event 
numbers of 10 x r, 18.5 x r and 21 x r at m a = 8 GeV, Mt 1s and 10.5 GeV, respectively. 
The statistical significances of the a peaks for L = 10 pb -1 are then rx the r = 1 results 
of 0.16(7, 0.08cr and 0.22(7, respectively. 

Of course, we currently expect that substantial early running will mostly take place at 
y/s = 7 TeV and ^fs = 10 TeV. As noted earlier, lower ^/s implies a somewhat smaller a 
cross section in the [8 GeV, 2m#] mass interval on which we are focusing. Roughly, relative 
to yfs = 14 TeV, the a cross section decreases by a factor of ~ 1.3 at yfs = 10 TeV and 
a factor of ~ 1.7 at y/s = 7 TeV in this mass interval. Since the backgrounds are also 
basically gg fusion induced, we presume that these same factors will apply to them. At 
^/s = 10 TeV (-^/s = 7 TeV) this then will reduce the statistical significances given above 
by a factor of l/y/1.3 (l/y/1.7). The statistical significances at m a = 8 GeV, My is and 
10.5 GeV are, respectively, then 0.14o", 0.07a, 0.19a at 10 TeV and 0.12a, 0.06cr, 0.17a at 
7 TeV, all to be multiplied by r. 

Given the above results, we can tabulate the integrated luminosity L needed to achieve 
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Table 2: Luminosities ( fb 1 ) needed for 5cr if tan/3 = 10 and cos 6a = 0.1. 



Case 


m a = 8 GeV 


m a = M Tl5 


m a <; 2m B 


ATLAS LHC7 


17/r 2 


63/r 2 


9/r 2 


ATLAS LHC10 


13/r 2 


48/r 2 


7/r 2 


ATLAS LHC14 


10/r 2 


37/r 2 


5.4/r 2 



a 5<t significance at each of the three energies. The results appear in Table ||. The required 
L's away from the Upsilon resonance may be achieved after a year or two of LHC operation. 
The sensitivity of the required luminosities to r shows the importance of firmly establishing 
the precise efficiencies for background and signal. We look forward to continued and 
detailed work by the ATLAS collaboration in this area. Of course, we must not forget 
that the required L's are very sensitive to tan/3, cos 9a and BR{a — > very roughly 

for tan/3 / 10, cos 9 a / 0.1 and/or BR(a — > ^ 0.003 the tabulated luminosities 

need to be multiplied by 

/ 0.003 \ 2 / 0.1 \ 4 / _10_\ 3 - 2 " 3 - 6 

\BR(a fi+fJ,-)) \cos6a) \tan/3y ' ^ ' ' 

where the 3.2 applies for m a ~ 8 GeV and the 3.6 applies for m a < 2m#. Depending upon 
the precise value of m a and tan/3, in the m a mass range under discussion Fig. p] shows that 
BR(a — > can range from a low of 0.0023 at tan/3 = 1.5 and m a ^ 2m# to a high of 

0.0033 for tan/3 > 3 and m a = 8 GeV. The minimum values of cos^, with and without 
placing a maximum on the light-a finetuning measure G, were detailed in Table |l|. 

Studies by CMS analogous to the ATLAS studies discussed above are under way |35| . 

5. Conclusions 

In this paper we have shown that a dedicated analysis of the dimuon spectrum at the Teva- 
tron and LHC at low masses, i.e. M^+^- < 2ms, will provide very important constraints on 
models containing a light CP-odd Higgs boson. We employed the published L = 630 pb _1 
CDF analysis of the dimuon spectrum between ~ 6.3 GeV and 9 GeV by CDF and found 
that constraints on the bba coupling C ab i become competitive with those from T(nS) — > 7a 
decays for 8.5 GeV <^ m a ^ 9 GeV, and will be superior for larger data sets. In addition, 
only hadron colliders have the kinematic reach to constrain |C a ^| in the important region 
Mt 3S ^ m a ^ 2mB- In particular, for L = 10 fb~ , the Tevatron will provide significant 
constraints on the |C ab ^| > 1 portion of the 8 GeV < m a < 2ms mass region that would 
allow an NMSSM ideal Higgs scenario with an ~ 100 — 105 GeV CP-even h decaying 
primarily via /i->aa-f t + t~t + t~ to be possible with neither electroweak finetuning nor 
"light-a" finetuning. It is also very noteworthy that our rough estimates of the limits that 
CDF could place on |C afe ^| using the L = 630 fb -1 event rates in the 9 < m a < 12 GeV 
region are such that the observed a M discrepancy could not be explained by a light a. 
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For the LHC, we have obtained rough estimates of what will be possible using infor- 
mation available from the ATLAS collaboration, in particular regarding the efficiency (for 
triggering, tracking, pt cuts, etc.) for retaining a — > events. We find that it will be 

possible to obtain a 5a signal for a light a with tan/3 ~ 10 and cos^ ~ 0.1 throughout 
the entire range 8 GeV <^ m a <~ 2m# away from the Upsilon peaks for L ~ 13 fb _1 at 

= 10 TeV or L ~ 17 fb _1 at -/s = 7 TeV. For example, at m a = 10.5 GeV, only 
L ~ 7 fb _1 at -v/i = 10 TeV or L ~ 9 fb _1 at ^ = 7 TeV is required to achieve a 5a signal 
for such an a. 

Of course, not all acceptable NMSSM models have {C^l as large as ~ 1. As an extreme 
example, at tan/3 = 1.7, cos^ ~ 0.1 is possible for small light-a finetuning (corresponding 
to CVit ~ 0.17). In this case, the a cross section at m a <J 2ms is about a factor of 18 
smaller than at tan/3 = 10 and cos6a ~ 0.1. Using statistical extrapolation this suggests 
that as much as 324 times more luminosity would be needed to achieve the same statistical 
significances as above. However, one should keep in mind that it may in the end be possible 
to obtain net efficiencies for the a at ATLAS and CMS in excess of the current ATLAS 
estimate of 10%. Indeed, early CMS studies suggest that net efficiencies might be as high 
as 30% . Since the needed L scales inversely with the square of the efficiency, assuming 
yfs = 14 TeV and r = 3 one finds that a 5a signal could be achieved for tan/3 = 1.7 and 
cos 6a ~ 0.1 with L ~ 195 fb _1 , an integrated luminosity that should be achieved in the not 
too distant future, although background levels might be larger at the higher instantaneous 
luminosities needed to achieve such large total L. 

Overall, this kind of search is quite important given that there are many models in 
which light a's are present that have significant, even if not enhanced, couplings to gluons 
via quark loops and that would have reasonable a —> ji + ^ branching ratio. Searching for 
such a's and constraining their possible masses and couplings is an important general goal 
and both Tevatron and LHC data will be of great value. 
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